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Precise comparisons of bottom-pressure and altimetric ocean tides 
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[i] A new set of pelagic tide determinations is constructed from seafloor pressure 
measurements obtained at 151 sites in the deep ocean. To maximize precision of estimated 
tides, only stations with long time series are used; median time series length is 567 days. 

Geographical coverage is considerably improved by use of the international tsunami 
network, but coverage in the Indian Ocean and South Pacific is still weak. As a tool for 
assessing global ocean tide models, the data set is considerably more reliable than older 
data sets : the root-mean-square difference with a recent altimetric tide model is 
approximately 5 mm for the M 2 constituent. Precision is sufficiently high to allow 
secondary effects in altimetric and bottom-pressure tide differences to be studied. The 
atmospheric tide in bottom pressure is clearly detected at the Si, S 2 , and T 2 frequencies. The 
altimetric tide model is improved if satellite altimetry is corrected for crustal loading by the 
atmospheric tide. Models of the solid body tide can also be constrained. The free core- 
nutation effect in the Ky Love number is easily detected, but the overall estimates are not as 
accurate as a recent determination with very long baseline interferometry. 

Citation: Ray, R. D. (2013), Precise comparisons of bottom-pressure and altimetric ocean tides, J. Geophys. Res. Oceans , 118, 4570- 
4584, doi: 10.1002/jgrc.20336. 


1. Introduction 

[ 2 ] Aiming for measurements of finer and finer precision 
is generally a prescription for progress in the natural scien- 
ces. For tidal phenomena, which result from forces across 
the entire Earth system, improved precision can reveal 
subtle, sometimes unexpected, effects. A now-classic 
example is precisely measuring tidally forced nutation and 
thereby discovering the excess flattening of the core-mantle 
boundary [Gwinn et al., 1986], This paper attempts to 
obtain improved complementary measurements of deep- 
ocean tides with satellite altimetry and with bottom- 
pressure recorders (BPR). The result is concomitant insight 
into atmospheric and solid-earth tides. 

[ 3 ] The original motivation for this work, it may be con- 
fessed, was to obtain better validation data for testing new 
global tide models, and thereby to update older work by 
Shum et al. [1997] and others. The data so obtained are 
indeed useful to that purpose ; such model comparisons will 
be addressed elsewhere within the context of a much larger 
effort with a multitude of different tests (D. Stammer, per- 
sonal communication, 2012). The present paper describes 
the construction and compilation of the improved pelagic 
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test data set and then explores what such improved preci- 
sion might reveal. 

[ 4 ] One of the standard and most widely used data sets 
for testing deep-ocean tide models was compiled in the 
early 1990s by David Cartwright and Christian Le Provost 
and their colleagues. That data set comprised tidal har- 
monic constants from 102 stations, 50 from open-ocean 
island tide gauges, and 52 from seafloor pressure recorders. 
The data set was used by Shum et al. [1997, Table 3] and 
many others since. Section 2 describes construction of a 
new “ground truth” data set, based strictly on bottom- 
pressure data, which is found to be considerably more accu- 
rate than the old 102-station set. The remaining sections 
address effects that can be seen in the altimeter minus 
bottom-pressure differences, with section 4 devoted to 
atmospheric tides and section 5 to solid-earth tides. 

[ 5 ] The paper addresses only tidal constituents of diurnal 
and higher frequency. No attempt is made to construct a 
reliable data set of long-period tides, even though there is a 
great need for one. In light of their longer periods and the 
higher background noise levels at those periods, a reliable 
compilation of station tide estimates would generally 
require considerably longer time series than those adopted 
here. Probably the most reliable test data set for long- 
period tides is still that compiled by D. Luther and pub- 
lished as an appendix in Miller et al. [1993]; those were 
based on multidecade time series from small Pacific Ocean 
islands. 


2. New Deep-Ocean “Ground Truth” Data 

[6] The new test data set consists of tidal harmonic con- 
stants at 1 5 1 stations, distributed about the globe as shown 
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Figure 1. Locations of the 151 deep-ocean bottom-pressure stations. 


in Figure 1. The stations are intentionally restricted to the 
deep open ocean, the traditional domain of satellite altime- 
try and where altimeter-based tide models are most accu- 
rate. Complementary test data sets restricted to shelf and 
near-coastal tides have also been developed [Lyard et al., 
2006; Ray et al., 2011] and continue to be expanded and 
improved, but these will not be used in the present work. 

[7] Tidal constants for a majority (86) of the 151 stations 
have been specially computed for this work from the meas- 
ured bottom-pressure time series (sampled hourly or 
faster). In addition, new tidal analyses have been computed 
for 14 time series by Doug Luther and Martin Guiles and 
for six time series by Chris Hughes, and these are also 
reported here for the first time with their permissions. 

2.1. Station Selection Criteria 

[8] The 151 stations were selected based on several cri- 
teria. Foremost, only tidal constants derived from seafloor 
pressure records are used. Bottom-pressure data offer a 
number of advantages over island tide-gauge data: (1) 
Within the tidal frequency band and above, the background 


noise in bottom pressure is almost always markedly lower 
than surface height measurements; Figure 2, comparing the 
sea-level spectrum at the island of Nantucket with bottom- 
pressure spectrum from a site about 200 km distant, is a 
very typical example. (2) Bottom pressure is less sensitive 
to internal tides, which is desirable, assuming it is a baro- 
tropic model being tested and not a fully three-dimensional 
model. (3) Island gauges can be influenced by more local- 
ized effects, such as pronounced amplitude or phase distor- 
tions not representative of the deeper surrounding ocean 
[e.g., Farrow and Brander, 1971 ; Douillet, 1998], 

[9] Second, station data were selected according to the 
length of pressure time series, with short time series 
rejected. Longer time series generally result, of course, in 
improved estimation of harmonic constants, and they also 
allow better separation of constituents with closely neigh- 
boring frequencies, such as the P 1 -S 1 -K 1 triplet which are 
separated in frequency by only 1 cpy. Figure 3 summarizes 
in histogram form the time series lengths for all selected 
stations. No selected series is shorter than 90 days, 80% of 
the series are of at least 1 year duration and 35% are 2 years 
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Figure 2. Spectrum of (a) sea level at Nantucket, USA, and (b) bottom pressure at DART station 
44402 (near 39.5°N, 70.6°W). The two stations are separated by roughly 200 km. Both spectra are based 
on approximately 1.5 years of data. While diurnal and semidiurnal tidal peaks are comparable at both 
sites, the bottom-pressure background energy in Figure 2b is much lower across the entire frequency 
range. 
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Figure 3. Histogram of lengths (in days) of all selected 
bottom-pressure time series. The shortest time series is 90 
days. The median is 567 days. 


or longer. It is only within the last few years that many of 
these very long time series have become available. Spencer 
and Vassie [1997] describe the first multiyear BPR time se- 
ries (station MYRTLE) from the deep ocean that did not 
entail repeated redeployments. 

[ 10 ] Third, stations located in close vicinity have been 
culled to avoid overweighting any small ocean area. This is 
only partially successful, and it is clear from Figure 1 that, 
for example, the South Atlantic is sampled far more heavily 
than the South Pacific. To some extent, the spacing of 
Topex/Poseidon altimeter tracks guided this culling of sta- 
tions, since for some of the estimations in sections 4 and 5 
the important criterion is that the sampling is independent. 
For example, Figure 4 shows the selected stations in the vi- 
cinity of the Drake Passage, and even though these stations 
are fairly close on a global or basin scale, the stations are 
still seen to sample the altimetry in a reasonably independ- 
ent manner, based on which tracks are closest to which 
stations. 

[n] This last criterion is the most subjective, but it is 
critical to avoid overweighting since some modem experi- 
ments deploy very tightly packed arrays. The KESS array 
[ Donohue et al., 2010], for example, could alone provide 
44 stations and would thereby overwhelm any global statis- 
tics. In cases like this, an attempt was made to maximize 
areal coverage by choosing a handful of the most widely 
separated stations. In many cases the choice of which 
nearby stations to select can be made based on time series 
length or some other judgment about data quality. A few 
examples will be noted below. 

2.2. Data Sources 

[ 12 ] Tidal harmonic constants derived from 414 bottom- 
pressure time series were previously collected in three pub- 
lications sponsored by the International Association for the 
Physical Sciences of the Oceans (1APSO): Cartwright et 
al. [1979]; Cartwright and Zetler [1985]; Smithson 
[1992]; along with additional data subsequently included 
in online archives at the U.K. National Oceanography 
Centre, Liverpool (formerly the Proudman Oceanography 
Laboratory (POL)). Many of these IAPSO time series, 
however, are either very short or represent redeployments 
at identical locations. For example, removing all series 


shorter than 200 days and all stations located within 100 
km of another station, as well as all stations in water shal- 
lower than 250 m, reduces the 414 time series to only 53, 
and many of these are near other, more modern stations. 
The IAPSO compilations are the source for 25 of our 151 
stations. The shortest time series evident in Figure 3 were 
all IAPSO stations, but they were retained to help improve 
geographical coverage. 

[ 13 ] An additional 15 stations, although part of the 
IAPSO compilation, have instead been extracted from the 
Global Undersea Pressure (GLOUP) archives, housed at 
POL. The GLOUP archive includes more extensive tidal 
analyses than those recorded in the older IAPSO booklets, 
which were generally restricted to eight major constituents, 
so the GLOUP data are therefore preferred for present pur- 
poses. Five GLOUP stations not included in the IAPSO 
compilations have also been added; they include two time 
series of length 730 days. GLOUP pressure data corre- 
sponding to IAPSO stations 1.1.68 and 1.1.69, being two 
successive occupations on the Azores Rise but separated by 
only 10 mi, have here been combined to form a 1 1 month 
time series and tidally analyzed. The GLOUP data include 
stations from the well-known ACCLAIM program deploy- 
ments in the South Atlantic and Southern Oceans [Spencer 
et al., 1993]. 

[ 14 ] Most of the IAPSO data, including many old analy- 
ses from the earliest days of seafloor records [e.g., Eyries et 
al., 1964; Snodgrass, 1968], are too short to meet our 
adopted criterion of 90 days duration. Fortunately, over the 
past two decades there have been a number of oceano- 
graphic campaigns collecting long series of pressure data, 
which can considerably augment the IAPSO compilations. 

[ 15 ] Major recent BPR campaigns worth explicitly men- 
tioning in this regard include : 

[16] 1. MOVE: an array of three stations near 16°N 
between the Lesser Antilles Arc and the Mid-Atlantic 


o 



Figure 4. Locations of selected bottom-pressure stations 
near the Drake Passage, with ground tracks of the Topex/ 
Poseidon and Jason satellites. The two westernmost stations 
are from unpublished data obtained courtesy of R. Watts, 
T. Chereskin, and K. Donohue. 
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Ridge [Kanzow et al., 2006] plus another station near 
21°N. The time series ranged from 591 to 809 days. 

[n] 2. ASTTEX: five stations along a Topex/Poseidon 
(T/P) satellite groundtrack, southwest of Namibia in the 
South Atlantic [Ray and Byrne, 2010]. Data for 11 stations 
are available along this track; every other station was 
selected here, with an average separation of 255 km. All se- 
ries are of about 800 days duration. 

[is] 3. AWI-ACC: seven selected stations in the South 
Atlantic, several of which extend the T/P line begun by 
ASTTEX, collected by the Alfred Wegener Institute 
(AWI). Six of these time series exceed 1000 days duration. 

[ 19 ] 4. KESS: seven stations selected from a dense array 
of 44 stations sitting in or near the Kuroshio Extension 
[ Donohue et al., 2010]. The seven were chosen to maximize 
areal coverage; the closest pair is separated by 168 km. 

[ 20 ] 5. HOME: nine stations surrounding the Hawaiian 
Ridge, collected by Alan Chave’s group at Woods Hole 
Oceanographic Institution between April 2001 and May 
2002 (about 380 days) and tidally analyzed by Doug 
Luther. 

[ 21 ] 6. RAPID: a total of nine selected stations from a 
wide-ranging array of stations in the North Atlantic, includ- 
ing a mid-ocean line at 24.5°N and several arrays off the 
U.S. east coast, deployed to study rapid changes in the 
Atlantic overturning [Rayner et al., 2011; Elipot et al., 
2013]. From six stations clustered very close to the African 
coast, just south of the Canary Islands, the two stations far- 
thest from land (EB3 and EBH2) were selected, with sepa- 
ration distance 162 km. Tidal constants for four of the 
selected RAPID stations were provided by Chris Hughes. 
He also provided constants for two stations in the 
Argentine Basin of the South Atlantic [ Hughes et al., 2007]. 

[ 22 ] 7. UR1-GOM: three stations from an array of sta- 
tions in the central [ Donohue et al., 2006] and western 
[Donohue et al., 2008] Gulf of Mexico. The three selected 
are the most widely separated of the time series that 
exceeded 1 year. Watts’s group (University of Rhode 
Island) also provided data from several other campaigns 
back to 1996. 

[ 23 ] 8. DART : a widely distributed set of 47 stations 
selected from the international (primarily American) tsu- 
nami early warning network, with 15 min data publicly 
available via the National Data Buoy Center. The DART 
data, originally recorded in units of pounds per square inch 
absolute (PSIA), are distributed in units of meters of equiv- 
alent water height after applying a uniform scaling of 670 
mm/PSIA conversion; for the present work the units have 
been converted back to pressure in Pa, taking 1 
PSIA = 6894.757 Pa. 

[ 24 ] One station — Cobb Seamount, taken from Larsen 
and Irish [1975] — is worth mentioning, because it is the 
one site violating the deep water criterion, but only because 
the recorder sat at the peak of a barely submerged seamount 
(depth 35 m) in an otherwise open expanse of the deep 
ocean. 

[ 25 ] Particularly important for the present work is the 
DART array, comprising nearly one third of the whole 
BPR data set. Its wide distribution of stations considerably 
strengthens coverage of the global ocean, especially fol- 
lowing the DART program’s extensive augmentation after 
the great 2004 Sumatra earthquake and tsunami. The new 



Figure 5. Monthly estimates of the amplitude and phase 
lag of the M 2 tide for DART station 46412, located west of 
San Diego, California. Vertical dashed lines mark times of 
redeployment of the bottom-pressure recorder. Data 
obtained in 2010 and afterwards were about 27 km distant 
from earlier data. The phase jump in late 2006 probably 
results from a time-tag error since a similar offset, but half 
as large, occurs for the Oi constituent. 

stations add valuable data in some previously unsampled 
regions. Use of the DART data, however, does require 
some care for at least two reasons: (1) There are a fair 
number of outliers, or generally suspect data, which require 
special handling. (These likely stem from inevitable prob- 
lems associated with real-time operations from the sea- 
floor.) Andreas Macrander (AWI), who has processed 
DART and other BPR data as part of validation exercises 
for satellite gravity missions [Macrander et al., 2010], has 
very kindly made available his corrected time series. Of the 
47 DART tidal analyses reported here, 36 are based on 
Macrander’ s time-series data. (2) Like any in situ instru- 
mentation in harsh environments, DART pressure recorders 
are periodically serviced and redeployed. The redeploy- 
ments are generally close to previous deployments for a 
given station, but sometimes insufficiently close to sample 
precisely the same tide. Figure 5 displays a possible exam- 
ple of this problem for DART station 46412, located off the 
coast of California. Shown is a series of monthly estimates 
of the M 2 tide for this station, derived by a response analy- 
sis, with vertical lines marking times of redeployments. 
According to available documentation, the May 2009 move 
entailed the largest change in position, about 27 km north- 
east toward the California coast, and the estimated ampli- 
tudes do seem to reflect this, although the slow rise in mid- 
2009 is curious and unexpected (it is unlikely to be from 
initial deployment stresses, since the instrument had 
actually been deployed in July 2008, overlapping the earlier 
deployment). The large phase jump in October 2006 corre- 
sponds to an insignificant change in horizontal position, 
only about 4 km. In fact, the 2006 phase jump of approxi- 
mately 7° is almost surely the result of a time-tag error of 
approximately 15 min; the phase of Oi jumped by a corre- 
sponding amount of 3 4’. For this DART station the 
adopted tidal constants are those computed for the period 
November 2006 through May 2009 (and obviously 
assigned to the old location). Station 46412 is not the only 
DART station with such problems. 

[ 26 ] A few DART stations were unused if they fell too 
near other stations. For example, DART station 51407 is 
close to Hawaii and unneeded in light of the HOME array. 
Station 55013, only 44 km from 55015 but with a shorter 
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time series, is similarly unneeded. On the other hand, 
44402 is preferred over RAPID station RS3 (separated by 
75 km) because it has a longer time series. 

[ 27 ] A small number of stations that passed initial accep- 
tance criteria were still rejected because the data were sub- 
sequently deemed suspect. For example, station M2. 5 of 
the MOVE array appears to have a timing error of roughly 
20 min, detected in part by comparison to nearby stations 
M2 and M3. But the station is unneeded in any event 
because of these nearby stations, so no loss of coverage 
arises by rejecting M2. 5. The most conspicuous outliers 
were the five stations in the 1986-1987 BEMPEX array of 
the North Pacific, which were based on constants reported 
by Filloux et al. [1991] and subsequently included in the 
1APSO catalog [Smithson, 1992]. The suspicious data have 
recently been traced to a flawed tidal analysis, and Doug 
Luther, who has the original time series, has kindly pro- 
vided constants from a new tidal analysis that corrects the 
earlier errors. 

2.3. Tidal Analysis of BPR Data 

[ 28 ] A total of 86 BPR time series were tidally analyzed 
as part of the present study. Hannonic constants for 25 sta- 
tions, either GLOUP stations or some RAPID stations, 
were computed by C. Hughes and others at the Proudman 
Laboratory. Constants for an additional 14 stations were 
obtained from Luther. The remaining 26 stations were from 
the IAPSO compilations or other published sources [e.g., 
Larsen and Irish, 1975]. 

[ 29 ] All stations that entail time series shorter than 160 
days were extracted from the IAPSO compilations. The 
tidal analyses of those time series were obviously per- 
formed by a variety of different investigators using a vari- 
ety of methods. The shortest time series were, very 
justifiably, analyzed by a response analysis, using either a 
reference tidal prediction, if available and appropriate, or 
the equilibrium tide [ Cartwright et al., 1969]. 

[ 30 ] All other time series were analyzed by least-squares 
harmonic methods, with standard adjustments for nodal 
and perigee variations within constituents. None of the 
time series was long enough, however, to support separat- 
ing third-degree from second-degree tides. Third-degree 
tides are generally reduced in amplitude relative to second- 
degree tides by a factor 1/60 — the lunar parallax — and 
rarely exceed 1 cm even in resonance [Cartwright, 1975], 
except in isolated gulfs that do not apply here (e.g., the 
third-degree N 2 constituent reaches 2.5 cm on the coast of 
Maine). Thus, the effects of inseparation should be minor 
and were deemed an acceptable tradeoff of using BPR data. 

[ 31 ] The number of constituents included in any station 
analysis varied according to both the length of time series 
and the estimation noise levels. In the distributed data 
accompanying this paper, a reduced set of common constit- 
uents is used for all tides, although not all are available for 
all stations. For example, of the 151 stations, only 121 
include constituent 2Qi. Selected compound tides have 
been included up to species 4. 

[ 32 ] Following each tidal analysis the estimated ampli- 
tudes were adjusted to account for any temporal averaging 
of the raw pressure measurements. This adjustment, often 
overlooked but necessary for high-precision work, arises 
because many of the BPR time series consist not of spot 


values of the pressure at discrete times but rather averages 
over some interval, often 1 h. Such averaging will slightly 
reduce the tidal amplitude of a constituent with period T by 
the amount (e.g., Malin and Chapman [1970], Appendix A) 

P/2 

— f cos(2nt/T)dt = (T/nP)sm(nP/T) 

-P/2 

where P is the averaging interval. Thus, when using hourly 
averages, semidiurnal amplitudes should be increased by a 
factor 1.0115, quarter-diurnal amplitudes by 1.0472. Tidal 
constants for the ASTTEX, KESS, the Gulf of Mexico, and 
older arrays from Watts’s group (e.g., from the Labrador 
Sea) have been so adjusted. Recent data from POL are gen- 
erally 1 5 min averages, which results in very little averag- 
ing error, so no adjustments have been done to these data. 
The DART data generally consist of 15 min spot readings, 
so no adjustment is warranted. The BEMPEX and HOME 
data were also analyzed (by Luther and Guiles) at 15 min 
sampling. Owing to lack of documentation, no corrections 
can be made to the old IAPSO data. 

[ 33 ] A second adjustment was done to attempt to remove 
from the bottom pressures the effects of atmospheric tides. 
This is discussed in more detail in section 4 below. To the 
extent that this removal is successful, the final bottom- 
pressure tidal constants then represent purely the ocean 
tide. 

2.4. Conversion of Bottom Pressure to Equivalent 
Surface Elevation 

[ 34 ] How tidal bottom-pressure fluctuations are to be 
converted to equivalent sea-surface heights for comparison 
with altimetry is a subject of some importance to the re- 
mainder of this paper. To the extent the subject is men- 
tioned in the tidal literature, the conversion is generally 
performed by dividing the pressure by gp mean , where p mean 
is the mean density of the water column and g is the local 
acceleration of gravity, and more often by simply employ- 
ing the conversion factor 1.01 mbar cm -1 [e.g., Smithson, 
1992], 

[ 35 ] A simple thought experiment reveals the inadequacy 
of this conversion. Consider a homogeneous ocean where 
density is solely a function of compression from the weight 
of the overlying water column. Then a tidal fluctuation at 
the surface of height r/ results in a water column identically 
compressed as before, save for a slab of water at the sea- 
floor of thickness p. The pressure fluctuation at the bottom 
thus corresponds to r/gp bottom . 

[ 36 ] The subject becomes more complex for an inhomo- 
geneous ocean in which internal tides occur. Since all tidal 
analyses here arise from long time series representing in 
essence temporal means of harmonic constants, only tem- 
porally coherent internal tides are potentially of impor- 
tance. At this stage, however, global maps of coherent 
internal tides — even for mode-1 M 2 — are too unreliable to 
allow us to account for the effects in any systematic way; 
the topic is an ongoing research concern [e.g., Dushaw et 
al., 2011] and, in fact, will probably require additional new 
data from future satellite missions [e.g., Fu et al., 2012]. 
Thus, internal tides are ignored in what follows. At some 
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Figure 6. Scale factor for converting bottom pressure to equivalent sea-surface height, in units of cm 
hPa -1 . (a) The standard conversion factor 1 /gp mean where p mean is the column-mean seawater density, (b) 
The conversion factor from equation (1). Note that in shallow seas the conversion factor approaches unity. 


future time, the influence from internal tides should be 
reconsidered, because it could well result in further 
improved precisions between BP and altimetric tides. 

[ 37 ] The following two paragraphs follow a line of 
thought suggested by C. Garrett (personal communication, 
2009). Consider a point in the ocean of depth H. The bot- 
tom pressure is gHp msm . Now consider a tidal conver- 
gence at this location, which acts to increase mass at 
every level of the water column by a small fraction e. The 
sea-surface height fluctuation rj would be eH if additional 
compression of the water column is ignored, but with 
compression it must be slightly less than this. The pressure 
p at every level, including the bottom, does increase by 
the same fraction e, the local density thus increases pro- 
portionally by the fraction epdp/dp = ep/c 2 , where c is 
the speed of sound. 

[ 38 ] The thickness of each level (i.e., distance between 
isopycnals) changes in two ways — increasing by e from the 
additive mass and decreasing by compression from the 
extra overlying mass; that is, the thickness changes by 
e(l — p/ pc 2 ). The second correction term is roughly gz/c 2 
at depth z and its average over the water column is gH /2c 2 , 
which provides a compression correction to the standard 
pressure-to-elevation conversion. Thus, for a bottom- 


pressure tidal fluctuation 6p h the corresponding surface ele- 
vation is 

V = s Pb{\ ~ gff/2c 2 )/fep mean ) (1) 

[ 39 ] The factor (1 — gH /2c 1 ) in deep water amounts to 
about a 1% adjustment, and less in shallower water, so it is 
small but still possibly important. The overall conversion 
factor rj/Spb that converts bottom pressure to equivalent 
sea-surface height is shown in Figure 6, where the relevant 
seawater densities (and the implied speed of sound) are 
extracted from the World Ocean Atlas and g is taken as 
normal gravity, dependent on latitude. The upper panel dis- 
plays the traditional factor l/gp mean and the lower panel 
displays the right-hand factor from equation (1). The larg- 
est differences arise in the deepest parts of the ocean where 
compressibility comes more into play. Note that in shallow 
seas the conversion factor approaches unity, 1 cm/mbar. 

[ 40 ] In fact, the altimeter-minus-BP tide differences are 
sufficiently precise to be sensitive to this method of con- 
verting between bottom pressure and equivalent sea- 
surface height. Root-mean-square (rms) differences 
between BPR tides and the GOT4.8 model, discussed in 
much more detail in the next section, are shown in Table 1 
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Table 1. RMS Differences (cm) Between BP Stations and Model 
GOT4.8 According to How Bottom Pressure Is Converted to 
Equivalent Elevation 3 


Conversion Method 

m 2 rms 

Oi RMS 

1 .00 mb/cm 

0.777 ±0.039 

0.411 ±0.020 

1.01 mb/cm 

0.586 ±0.034 

0.357 ±0.018 

SPmsan 

0.517 ±0.034 

0.323 ±0.019 

Equation (1) 

0.510 ±0.033 

0.296 ± 0.020 


“Uncertainties are standard deviations computed via bootstrap 
resampling. 


as a function of the method used to effect the pressure con- 
version. For both M 2 and Oi constituents the smallest rms 
values are obtained when using equation (1), with Oi show- 
ing clearest improvement over the more standard factor 
1 /gp m ean • The simpler conversions, such as a constant 0.99 
cm/mbar factor, result in considerably higher nns differen- 
ces, and they should probably be avoided except for use as 
a quick rule of thumb. 

[ 41 ] Finally for the special case of a homogeneous, but 
compressible, fluid that is neutrally stable, the potential 
density is constant so the density lapse rate is strictly the 
adiabatic lapse rate, dp/dz = —pg/c 2 . Thus, p(z) = p s ex p 
(— gz/c 2 ) ss p s (l — gz/c 2 ), where p s is the surface density 
and the approximation is valid because (gz/c 2 ) is generally 
of order <9(10~ 2 ) or smaller. Then 


Pmean = 1 + gH /2c 1 
Pbottom 1 f / C~ 


\-gH/2c 2 


so that (1) yields dpi, = pgPbottom > as noted above. 

3. Comparisons of BPR and Altimetric Tides 

[ 42 ] Table 2 is the main result of comparing the new 
BPR tide data set with the global altimetric tide solution 
GOT4.8. (Solution GOT4.8 and some variants of it are 


Table 3. RMS Differences (cm) With GOT4.8 



0, 

K, 

n 2 

m 2 

s 2 

Old 102-station data set 

0.77 

1.02 

0.64 

1.45 

0.83 

New 151-station data set 

0.30 

0.42 

0.25 

0.51 

0.37 


described in Appendix A. Since the altimetry is sensitive to 
the geocentric tide, the solid-earth contribution must be 
removed to yield a true ocean tide ; details for this are dis- 
cussed in Appendix B.) 

[ 43 ] The root-mean-square (rms) differences of Table 2 
were computed from all in-phase and quadrature compo- 
nent station-model differences, where model values were 
evaluated by bilinear interpolation at the locations of the 
151 BPR stations. These differences are tabulated for 18 
constituents. The table also includes a second statistic, the 
median of absolute differences, which is less sensitive to 
outliers (although such sensitivity can be useful for testing 
models or in other applications). 

[ 44 ] The quoted uncertainties in Table 2 represent an 
attempt to gauge statistical variability of the tabulated val- 
ues. Such assessment is necessary when, for example, com- 
paring different rms values from different models or 
different corrections. The uncertainties here are based on a 
bootstrap resampling of the station data (with replacement), 
which seems preferable to relying on dubious distributional 
assumptions. The tabulated quantities represent 1 standard 
deviation from these bootstrap resamplings. 

[ 45 ] The agreement between the BPR and altimetric tides 
must be considered quite good, well below 1 cm for all 
constituents. The largest rms of 5.1 mm coincides with the 
largest constituent, M 2 , while smaller constituents (identi- 
fied by the rms “signal” column in the table) show better 
agreements. 

[ 46 ] The quality of the new BPR data set is brought out 
clearly in Table 3, which compares nns differences (again 
with GOT4.8) for the five largest tides against similar cal- 
culations made with the old standard 102-station data set, 


Table 2. Comparisons Between Model GOT4.8 and Bottom-Pressure Tides 3 


Tide 

Number Stations 

BPR Signal RMS (cm) 

RMS Diff. (cm) 

Median Abs.Diff. (cm) 

2Q, 

121 

0.2 

0.085 ± 0.007 

0.040 ± 0.004 

Qi 

148 

1.8 

0.164 ±0.007 

0.107 ±0.008 

Oj 

151 

8.8 

0.296 ± 0.020 

0.181 ±0.012 

Pi 

146 

4.0 

0.234 ±0.011 

0.146 ±0.008 

Si 

117 

0.4 

0.331 ±0.029 

0.187 ±0.014 

K, 

151 

12.5 

0.423 ±0.016 

0.303 ±0.015 

h 

125 

0.7 

0.178 ±0.009 

0.1 10 ±0.008 

O 

O 

124 

0.4 

0.265 ±0.013 

0.151 ±0.013 

2N 2 

123 

0.8 

0.271 ±0.014 

0.146 ±0.014 

ft 

123 

1.0 

0.409 ± 0.028 

0.226 ±0.021 

n 2 

151 

6.4 

0.252 ±0.013 

0.156 ±0.007 

^2 

121 

1.1 

0.075 ± 0.003 

0.048 ± 0.003 

m 2 

151 

30.2 

0.510 ±0.033 

0.251 ±0.022 

L 2 

124 

0.8 

0.250 ± 0.020 

0.124 ±0.010 

T 2 

118 

0.6 

0.130 ±0.007 

0.075 ± 0.008 

s 2 

151 

11.2 

0.369 ± 0.022 

0.235 ±0.012 

K 2 

146 

3.1 

0.209 ±0.011 

0.129 ±0.011 

M4 

130 

0.2 

0.089 ± 0.004 

0.062 ± 0.004 


“Uncertainties are standard deviations computed via bootstrap resampling. 
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Difference (cm) 

Figure 7. Histogram of in-phase and quadrature compo- 
nent differences between bottom-pressure stations and the 
GOT4.8 model, for the M 2 constituent. The standard devia- 
tion a corresponds closely to the rms value quoted in Table 
2. Parameters 7 j and y 2 are skewness and excess kurtosis, 
respectively; the latter suggests the distribution is “long- 
tailed” with possible outliers. In fact, this is clear from the 
presence of several differences exceeding 3 < 7 . 


mentioned in section 1. The older data set has been used 
widely over the past two decades, especially for testing 
new tidal models, but Table 3 suggests it can now be con- 
sidered obsolete. Note in particular that the rms for M 2 
drops by a factor of 3 for the new test data set. As men- 
tioned above, about half of the old data set is based on 
island tide-gauge stations that may not properly reflect the 
open-ocean tide seen by altimetry. Island data may also be 
perturbed by internal tides, which are for the most part fil- 
tered out of altimetric tide solutions like GOT. Further- 
more, about one fifth of the old data set is based on time 
series of duration 1 month or less, and thus inherently 
much less precise than the tidal estimates compiled here. In 
light of all these reasons, it is perhaps not surprising, but 
nonetheless still satisfying, to see the rms for M 2 drop from 
1.5 cm to 0.5 cm. 

[ 47 ] A more detailed view of the station-model tide dif- 
ferences for M 2 can be seen in Figures 7 and 8. A histo- 
gram of the in-phase and quadrature component differences 
is shown in Figure 7. Relative to a normal distribution, the 
residuals are somewhat leptokurtic, suggestive of the pres- 
ence of outliers, which could be in either the station data or 
the altimeter estimates or both. It is relevant to note, how- 
ever, that the four stations with largest residuals are all 
from the old 1APSO compilation (specifically stations 
1.1.71, 1.2.107, 2.1.17, and 4.1.4). In fact, if the BPR sta- 
tions are partitioned into subsets, the M 2 rms difference for 
the 26 IAPSO stations is 0.79±0.11 cm, but only 

0. 34 ±0.04 cm for the 47 DART stations. 

[ 48 ] A standard scatterplot diagram displaying station 
versus model amplitudes or phases is not especially inform- 
ative because at fall scale nearly all points appear to fall on 
a straight line of unit slope. It is more revealing instead to 
plot amplitude or phase differences as a function of ampli- 
tude, which is shown in Figure 8 for M 2 and Figure 9 for 

01. While all differences are small, in keeping with the 
above statistics, Figure 8 reveals a slight tendency for 
GOT4.8 amplitudes to exceed BPR amplitudes above about 


45 cm. This discrepancy does not appear to be a simple 
scale error — for which many possible explanations might 
be offered — but more like a simple offset, which is harder 
to explain. The mean amplitude difference for all stations is 
0.77±0.41 mm, but for amplitudes greater than 45 cm the 
mean difference is 2. 8 ±0.8 mm. A similar discrepancy in 
large amplitudes is seen in N 2 (not shown). However, the 
tendency in 0[ appears opposite (Figure 9), so perhaps 
both are statistical flukes. The positive offset near very 
small (1 cm) amplitudes probably owes to additive noise in 
the GOT4.8 solution. There are no statistically significant 
mean differences in phase in either diagram; larger phase 
scatter for small amplitudes is, of course, expected and of 
no great concern. 

[ 49 ] Another curiosity can be noticed in Table 2: the 
model-station differences for // 2 are anomalously large. 
Now the GOT solutions for fj 2 , as well as the neighboring 
2N 2 , are not direct estimates but rather are extrapolated 
from the tidal admittances at N 2 and M 2 . Yet /x 2 seems rela- 
tively more error prone than 2N 2 even though the extrapo- 
lation is over a comparable frequency range (about 1 cycle/ 
month below N 2 ). In shallow water /x 2 is often significantly 
perturbed by M 2 -S 2 nonlinear interactions, which inhibits 
inference calculations, but the effect is expected to be slight 
in deep water. At Honolulu, however, Munk and Cart- 
wright [1966] noticed a substantial offset of the semidiur- 
nal admittance at /r 2 . Perhaps like M 4 in the deep Atlantic 
\Lyard et al., 2006; Ray, 2007], a significant part of /x 2 is 
nonlinearly generated in shallow waters and freely propa- 
gates into the deep ocean. If true, this suggests that relying 
on admittance calculations to map /x 2 may be too simplistic 
and that more direct estimation from altimetry and hydro- 
dynamics may be necessary. 

[ 50 ] In the diurnal band the GOT fields for 2Qi, Jj, and 
OO 1 are also all extrapolated from admittances, but their 
errors appear fairly small. The rms for OOi is larger than 


o o ^ o A, c*4b cO 0 

- gb-^oo O „ 0,0 o CO 


M2 



Amplitude (cm) 


Figure 8. Differences between BPR station tides and 
model GOT4.8 for M 2 (top) amplitudes and (bottom) 
phases as a function of station amplitude. Differences are 
in the sense: (model minus station). There are 151 differ- 
ence pairs in each panel. For amplitudes greater than 45 
cm, there is a tendency for GOT4.8 amplitudes to slightly 
exceed the BPR amplitudes. 
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Figure 9. As in Figure 8 except for tidal constituent Oi 


J i, even though its amplitude is generally smaller, but this 
is not unexpected since 00 1 is being extrapolated over a 
greater range of frequency. 

[ 51 ] Finally, the differences in Table 2 for Si also look 
large in relation to the constituent’s generally small size. 
This too is not unexpected. Si is notoriously difficult to 
measure — a multitude of systematic errors occur at the fre- 
quency of precisely 1 cycle/day — and because its proxi- 
mate cause is the Si atmospheric tide it is subject to a fair 
amount of temporal variability. 

[ 52 ] Aside from testing and comparing ocean-tide mod- 
els, the new bottom-pressure data set can be used to exam- 
ine interesting secondary features. The next two sections 
offer two examples, exploiting the fact that the BPR data 
are sensitive to atmospheric tides while altimetry is not 
(section 4) while the altimetry is sensitive to the solid-earth 
tide while the BPR data are not (section 5). 

4. Detection of Atmospheric Tides 

[ 53 ] As noted earlier, the BPR tidal data have been cor- 
rected to remove the effects of atmospheric tides. This sec- 
tion analyzes the extent to which this removal is both 
necessary and successful. Atmospheric tides impact the 
bottom pressure in two ways: (1) they directly add to the 
measured pressure by as much as 1-2 hPa and (2) they 
induce a dynamic ocean response, the so-called radiational 
tide. The latter is rightly considered part of the ocean tide 
and is not to be removed. Removal of the direct effect 


requires a model, either numerical or analytic, of the atmos- 
pheric tide. 

[ 54 ] The dominant barometric tide over the oceans is the 
solar S 2 . For this constituent, three models are here tested. 
The simplest is an analytic model given by Haurwitz and 
Cowley [1973], which is dependent only on latitude (j) and 
time : 


= l-161cosVsin(2f + 159 deg ), 

in hPa, where t' is local time in appropriate units. A sec- 
ond model is based on an analysis [Ray and Ponte, 2003] 
of 6 h surface pressure data in the operational analysis 
product of the European Centre for Medium-Range 
Weather Forecasts (ECMWF). Because 6 h sampling is 
inadequate to define S 2 , constraints were invoked concern- 
ing phase propagation to yield a realistic tide, although 
the result was possibly too zonally symmetric (see cited 
paper for details). Comparison with barometric station 
data showed reasonably good agreement, but it did reveal 
a 20 min timing discrepancy, which was corrected in the 
final phase data. A third model is also based on tidal anal- 
ysis of ECMWF operational surface pressure data, but 3 h 
rather than 6 h data, thus negating the need for assump- 
tions about phase propagation. These results were kindly 
provided by Jean-Paul Boy (University of Strasbourg). 
Analysis of both ECMWF models included annual modu- 
lations, for both Si and S 2 , which yields equivalent con- 
stituents Pi, Ki, R 2 , and T 2 . 

[ 55 ] Table 4 presents the main results of adjusting the 
BPR data for atmospheric tides, again using GOT4.8 as 
the comparison altimeter-based model. For these tests, a 
couple of points about the altimeter processing are rele- 
vant. Before tidal analysis the altimeter data were cor- 
rected for atmospheric loading, but only at frequencies 
below the tidal bands. This ensures that all tidal signals in 
the altimetry, including the radiational tides, are left intact 
and uncorrupted by any tidal signals in the air pressure 
fields. Secondly, some of the altimetric tidal constituents 
of Table 4 were not estimated directly but were inferred 
through admittance relationships based on the gravita- 
tional potential. For example, the T 2 constituent was 
inferred from S 2 under the assumption that radiational 
forcing is minimal in both, which to first order is probably 
satisfactory. Constituent R 2 cannot be so inferred, how- 
ever, since its gravitational forcing is so weak that radia- 
tional forcing is often dominant. Thus, R 2 is not included 
in the table. 


Table 4. Air-Tide Tests: RMS Differences (cm) for Model GOT4.8 a 


Air-Tide Model 

Pi 

Si 

Ki 

t 2 

S 2 

k 2 

S 2 (4.8a) 

None 

0.219 

0.416 

0.275 

0.157 

1.010 

0.209 

0.921 

Analytic 1 * 





0.372 


0.332 

ECMWF 6 h operational, 1 986-1 998 c 

0.241 

0.307 

0.280 

0.140 

0.376 

0.209 

0.334 

ECMWF 3 h operational, 2000-2006 d 

0.232 

0.303 

0.258 

0.139 

0.371 

0.208 

0.339 

Bootstrap a 

0.022 

0.022 

0.018 

0.015 

0.023 

0.016 

0.023 


“Based on 39 BPR stations in tropical latitudes. 
b From Haurwitz and Cowley [1973] 

“From Ray and Ponte [2003], 

d Air tides computed by J.-P. Boy, personal communication. 
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Table 5. Love Number h 2 (Real Component) 


Authors 

Method 

m 2 

0, 

K, 

Dehant et al. [1999] 

Theory 

0.6106 

0.6046 

0.5222 

IERS Conventions 2010 

Theory 

0.6078 

0.6028 

0.5236 

Herring and Dong [ 1 994] 

VLBI 

0.604 ± 0.002 



Ray et al. [1995] 

Altimetry 

0.613 ±0.007 

0.594 ±0.013 

0.512 ±0.012 

Haas and Schuh [1996] 

VLBI 

0.600 ±0.001 

0.606 ± 0.002 

0.496 ± 0.002 

Haas and Schuh [1997] 

VLBI 

0.600 ±0.001 

0.612 ±0.002 

0.512 ±0.003 

Wu et al. [2001] 

SLR 

0.606 ± ? 

0.618 ±? 

0.502 ± ? 

Rutkowska and Jagoda [20 1 0] 

SLR 

0.6151 ±0.0008 



Krasna et al. [2013] 

VLBI 

0.6072 ± 0.0003 

0.6026 ± 0.0009 

0.5267 ±0.0007 

This paper 3 

Altimetry 

0.611 ±0.003 

0.610 ±0.006 

0.540 ± 0.006 b 


“Listed uncertainties for this line correspond to one standard error. 
b After correcting altimetry for FCN effect in load tide; otherwise 0.543. 


[56] The rms values of Table 4 are based only on 39 sta- 
tions in tropical latitudes, since the air tide is very small 
outside this band. The air tide is very clearly detected in the 
altimeter-BPR differences for the main tides S 2 and S | . The 
rms reduction for S 2 is especially pronounced, dropping 
from 1.01 cm to 0.37 cm. The air tide is marginally 
detected in the annual sidelines T 2 and Kq (at least for the 
3-h ECMWF model). It is not detected at P b since the rms 
values increase for both ECMWF models. For all constitu- 
ents but one, and especially for K,, the lowest nns values 
are produced by the 3 h version of the ECMWF model. 
However, in light of the bootstrap-based standard errors, 
none of the differences between ECMWF models is signifi- 
cant, except for Kq for which the reduction in rms does just 
exceed 1-cr. 

[ 57 ] The final column of Table 4 uses a variant of 
GOT4.8, discussed in Appendix A, in which the altimetry 
is corrected for crustal loading by the S 2 atmospheric tide, 
a correction that has been hitherto always ignored. It 
shows a further reduction in the rms, suggesting the cor- 
rection is now warranted. See Appendix A for further 
discussions. 

[ 58 ] Finally, an attempt was made to detect the lunar 
barometric tide by using the same approach with a 
model based approximately on Haurwitz and Cowley 
[1969], The lunar tide, however, never exceeds 0.07 
hPa, and so it is probably unsurprising that the effect 
on altimeter-BP differences was negligible and the tide 
was not detectable. 


5. Estimation of Solid Body Tides 

[ 59 ] One possible contributor to differences between the 
altimetric tides and the bottom-pressure tides are errors in 
the adopted model of the solid body tide used to correct the 
satellite altimeter data. Thus, an analysis of altimeter-BP 
differences can conceivably be used to constrain models of 
the body tide and to estimate body-tide Love numbers. An 
early attempt to do this [Ray et al., 1995] was encouraging 
and found to be competitive with Love number estimates 
obtained from very long baseline interferometry (VLBI) 
[e.g., Mitrovica et al., 1994]. Our current BPR data set is 
much more accurate than that earlier tide-gauge data set. It 
appears, however, that VLBI data and analyses have 
improved considerably more. In principle, one expects 


VLBI to be the preferred method, since it measures the 
solid tide without the possible masking effects of inconsis- 
tent (and much larger) ocean tide signals in the altimeter- 
BP differences. 

[ 60 ] Let ( 9 , ip ) be spherical polar coordinates and t be 
time. The expression for the body tide of any constituent of 
degree-2 and species m (1 for diurnal waves, 2 for semi- 
diurnal), with frequency u>, is given by 

b(9, ip, t) = h 2 H cos(u!t + my))?™ (cosd) (2) 

where P 2 (cos9) is the associated Legendre function. The 
amplitude H is taken from the harmonic expansion of 
the tidal potential [Cartwright and Edden, 1973]; for 
example, for M 2 we have H = 8.1367 cm (taking 
63.194 cm from the Cartwright-Edden tables and scaling 
by their normalization factor (5/967r)^ 2 ). The Love 
number h 2 is taken as strictly real because (it turns out) 
standard errors are comparable or larger than anticipated 
imaginary components, which for M 2 is thought to be 
roughly 0.002-0.003, corresponding to an anelastic lag 
of 0.2° [Mathews et al., 1997]. 

[61] Table 5 presents a compilation of theoretical and 
empirical estimates of the real part of h 2 , including a value 
recommended by the International Earth Rotation and Ref- 
erence System Service (1ERS) Conventions [Petit and 
Luzum, 2010], the latter being mostly consistent with theo- 
retical calculations by Mathews et al. [1997]. There is a fair 
amount of scatter among estimates at the level of a few per- 
cent, yet all clearly display the free core-nutation effect 
that arises for Ki as its Love number is suppressed relative 
to nominal values [e.g., Wahr, 1981]. 

[ 62 ] Recent VLBI estimates of Krasna et al. [2013] 
agree closely with IERS values, lending strong evidence 
that these should now be considered definitive. The new 
altimeter-based estimates here appear not in very good 
agreement with these and must be considered much less 
accurate than the new VLBI estimates. In fact, discrep- 
ancies among all empirical determinations of h 2 , as seen 
in Table 5, suggest a long history of underestimating 
uncertainties, or at least underestimating systematic 
errors. It is therefore worth considering possible system- 
atic errors that may bias the altimeter-BP estimates. The 
main ones are probably errors in the load tide and in the 
satellite orbit. 
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[63] Consider first possible errors in the load tide. The 
component of the load tide of relevance here is the one 
induced by the degree-2, order- m prograde component of 
the ocean tide : 

C tm{9, <P, t ) = Dl m cos(ut + mip — ipl m ^P "2 (costf), (3) 

with amplitude Dj m and phase lag ip 2m . Dropping common 
factors of Pi, and time, we may write the combined 
in-phase body plus load tide, to which the altimeter-BP 
data are sensitive, as 

b + L 2fn = h 2 H + ah! 2 Dj m cosip^ (4) 

where a = (3/5)(p w /p e ) for mean densities of seawater 
and earth, p w and p e , respectively. The value of a is 
approximately 0.1125, which is probably accurate to ±1%. 
For M 2 the ocean-tide coefficients according to GOT4.8 are 

D 2m cos = 3.226cos 129.80deg = -2.065 ± 0.017cm 

where the uncertainty is taken from Egbert and Ray 
[2001]; uncertainties are somewhat larger for diurnal con- 
stituents. Our adopted It 2 from Farrell [1972] is —1.001; 
other published values for h 2 do not differ by more than 
2% \Nakada and Lambeck, 1987; Pagiatakis , 1990; Wang 
et al., 2012]. With nominal values of the Love numbers, 
the right-hand side of (4) for M 2 is then (4.95 cm + 0.23 
cm). So the load tide (for this degree-2 prograde compo- 
nent) amounts in total to about 5% of the body tide, and 
errors in it from errors in the ocean tide as well as errors in 
the various geophysical constants cannot affect body-tide 
estimates by more than about 0.1%. 

[64] An important exception, however, concerns K b 
where the loading number h 2 displays the same core reso- 
nance effect as seen in h 2 , but only for this one spherical 
harmonic component of the load. From Wahr and Sasao 
[1981] the perturbation in h 2 at the frequency of is 
about 15%, so failing to account for this perturbation can 
lead to errors in the estimation of h 2 of as much as 1%. In 
fact, accounting for this perturbation in diumal-band tidal 
loads is rarely done and was not done in GOT4.8. In light 
of this, the altimetric fields of used to produce the esti- 
mates of Table 5 have been recorrected for ocean tide load- 
ing to allow for this resonance; it reduces the size of the 
estimated h 2 although not nearly enough to bring it into 
alignment with the IERS Conventions. It is quite possible 
that other published estimates of Love numbers, including 
those using VLBI, also suffer from these errors in diumal- 
band load-tide corrections. 

[ 65 ] Another possible source of bias in the altimeter-BP 
Love numbers is satellite orbit error. The primary orbital 
ephemerides for T/P and Jason satellites were computed 
using tracking data from satellite laser ranging (SLR) and 
Doppler DORIS data with a heavy reliance on dynamical 
modeling \Lemoine et al., 2010]. Errors in the tide model 
adopted for the dynamical modeling can lead to tidally 
coherent errors in the orbit, and in the resulting sea-surface 
heights and thus in the estimated ocean tides, in a way dis- 
cussed by Colombo [1984] and Bettadpur and Panes 
[1994]. For the old Geosat ephemeris discussed by 


Bettadpur and Eanes the errors were quite large, reaching 
even 5 cm in some locations. Errors are certainly smaller, 
perhaps even by 2 orders of magnitude, for the modern tide 
models used for T/P and Jason, and we are here interested 
in only those errors that affect D 2 „,cos^J m , but the errors 
are surely nonzero. It is difficult to estimate this effect 
without knowing the true ocean-tide error. Unfortunately, 
estimated Love numbers are sensitive to errors of 1 mm or 
smaller in the altimetric-inferred body tide, so the altimeter 
estimates of h 2 may be affected. VLBI is of course immune 
to such errors ; the SLR estimates would not be. 

[66] Orbits based on GPS tracking are less sensitive to 
dynamical modeling errors, so an altimetric tidal solution 
built on such orbits (which could currently comprise Jason- 
2 and part of the Jason- 1 data) would be of interest. 

6. Summary Discussion 

[ 67 ] The new compilation of 151 precise estimates of 
bottom-pressure tides should prove a useful tool for assess- 
ing and comparing global ocean-tide models, and such 
work is in progress and will be reported elsewhere 
(D. Stammer, personal communication, 2012). The new 
data set is clearly superior to previous global compilations, 
as Table 3 emphasizes; rms differences for M 2 have 
dropped from 1.5 cm to 0.5 cm. This added precision 
should prove useful as models are further improved for the 
most stringent applications (which often involve correcting 
geodetic measurements of various types, including 
satellite-to-satellite tracking). 

[68] Because bottom-pressure recorders sense the atmos- 
pheric barometric tide, but satellite altimetry does not (at 
least not directly), the altimetric and bottom-pressure dif- 
ferences can detect this tide, very clearly for constituents 
S| and S 2 , and more marginally for T 2 and Kj. Our data 
may in future contribute to testing models of the atmos- 
pheric tide, since standard meteorological stations used for 
such testing are confined to only land [e.g., Ray, 2001]. 

[ 69 ] Similarly, since satellite altimetry directly senses 
the geocentric tide, altimeter-BP tide differences reveal the 
solid-earth tide, both the body tide of spherical harmonic 
degree 2 and the more spatially complex load tide. Deter- 
mination of the body-tide Love numbers is credible in the 
altimeter-BP differences for major constituents (Table 5), 
but they appear less precise than a recent VLBI determina- 
tion [Krasna et al., 2013], which agrees closely with the 
IERS Conventions for h 2 . The altimeter approach may be 
less accurate owing to some subtle systematic errors, dis- 
cussed in section 5, and also because VLBI is not affected 
by the additional large, and possibly inconsistent, ocean 
signals inherent in the altimetry and bottom-pressure data. 

[ 70 ] Nonetheless, detection of the body tide in altimeter- 
BP tide differences is more straightforwardly identified 
than the load tide, since the former depends on only one 
spherical harmonic component. No attempt has been made 
here to isolate a possible load-tide signal in the difference 
data. Such an attempt would benefit from using a dense 
array of stations, of the sort now being employed on land 
with dense arrays of GPS stations [e.g., Ito et al., 2009; 
Yuan and Chao, 2012], In fact, the full KESS array [Dono- 
hue et al., 2010], not decimated to a few stations as here, 
should be examined in this way, although, as with the 
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VLB1 estimates of Love numbers, we must expect that 
GPS data will be superior since GPS has no large ocean 
signal to contend with. 

[ 71 ] Aside from attempts to study air or solid tides, the 
precise altimeter-BP differences already point to certain sys- 
tematic errors that have previously been ignored in satellite 
altimetry. One is crustal loading by the S 2 atmospheric tide ; 
this signal has previously been absorbed into all ocean tide 
solutions and although small (1.5 mm) it is systematic and 
clearly detectable (Appendix A). Another systematic error 
involves inadequate modeling of the fluid core resonance 
effect in diumal-band loading Love numbers, the theory for 
which was laid out by Wahr and Sasao [1981], This error 
has generally been neglected in almost all geodetic applica- 
tions, but it has a not insignificant contribution to the present 
Love number estimate for K[ and should therefore be con- 
sidered in other studies requiring accurate ocean tide loading 
corrections, such as those using VLB1 or superconducting 
gravimetry. A third systematic error that is clearly detectable 
in the altimeter-BP differences concerns the presence of the 
atmospheric tide in the dry-tropospheric range correction; 
while this is allowed for in Jason GDRs, it was previously 
neglected in most T/P data processing. 

[ 72 ] It is somewhat astonishing to see such tiny 
millimeter-level signals come into play. But as noted at the 
beginning of this paper, it is often down in the limits of 
detection that interesting new phenomena are identified. 
The new precise tide estimates compiled here help push 
forward into those limits by confirming that some of these 
tiny signals are indeed real. 

Appendix A : GOT Altimetric Tide Solutions 

[ 73 ] The altimeter-based tidal constituent data sets 
employed in the main part of this paper are taken from the 
so-called Goddard/Grenoble Ocean Tide (GOT) series. 
(The “Grenoble” in the acronym stems from early versions 
that used a finite-element hydrodynamic model from Chris- 
tian Le Provost’s group — at that time located in Greno- 
ble — as a prior.) This appendix briefly describes the most 
recent GOT4 versions. 

[ 74 ] All of the GOT series are based on empirical tidal 
analyses of multisatellite altimeter data. The methodology 
essentially follows Schrama and Ray [1994] and Ray 
[1999], As in those works, tidal analyses were computed on 
altimetric residuals relative to an adopted prior model, which 
in the most recent versions was a combination of several 
global, regional, and local tide models blended across their 
mutual boundaries. In deep water standard hannonic analy- 
ses were used to estimate constituents Qi, Oi, Si, Ki, N 2 , 
M 2 , S 2 , K 2 , and M 4 . Smaller constituents (e.g., Ji in Table 2) 
were inferred from admittances. The computed amplitude 
and phase data are on a 0.5° global grid, which is more than 
adequate for the deep ocean but only marginally adequate in 
shelf seas and quite inadequate in coastal waters. 

[ 75 ] GOT4.7 has been used in a number of applications, 
so it is worth beginning discussion with that version. In the 
deep ocean between latitudes ±66°, which covers the do- 
main of primary interest here, GOT4.7 was purposely based 
on measurements solely from Topex/Poseidon (T/P) and no 
data from Jason- 1 or its follow-ons. The T/P and Jason sat- 
ellites, it turns out, are slightly inconsistent at the period of 


the S 2 tide, owing to some subtle and not completely under- 
stood effects dependent on solar radiation and the angle 
between the satellite orbit plane and the earth-sun axis, 
which varies at the alias period of S 2 (59 days). It is thus 
useful for the following tests to have a tide solution derived 
solely from T/P data. GOT4.7 did also use data from other 
altimetric satellites, including ICESat, but only in regions 
outside the domain of our new BPR test data set. The altim- 
eter data were adjusted for atmospheric loading by adopt- 
ing an inverted-barometer model, but only at periods 
longer than 2 days; this ensures that tidal signals are not 
corrupted in any fashion by this adjustment. 

[ 76 ] Some experimental GOT solutions have since been 
computed, which are nearly identical to 4.7, except for a 
few technical details. Most of these were meant to address 
problems with the solar constituent S 2 . In fact, GOT4.8 and 
GOT4.9 are identical to 4.7 except for S 2 (and its inferred 
neighbor T 2 ). 

[ 77 ] GOT4.8 corrects a problem with the dry- 
tropospheric correction that has been used in T/P Geophysi- 
cal Data Records (GDRs) since the beginning of the mis- 
sion. These T/P tropospheric corrections have been derived 
from 6-h ECMWF operational surface pressures. Because 6 
h sampling is inadequate to capture the full S 2 atmospheric 
tide signal, there is a corresponding error in the tropo- 
spheric correction at this frequency. The error is very small. 
The dry-tropospheric delay is given approximately by the 
surface pressure scaled by —2.277 mm/hPa [Hopfield, 
1971], and since the air tide never exceeds 1.2 hPa over the 
ocean, the tropospheric tide error never exceeds 3 mm. It 
is, however, systematic, of hemispheric scale. For GOT4.8, 
this tide error in the dry-tropospheric correction was mostly 
removed by using the approach advocated by Ponte and 
Ray [2002], which requires computing 6 h climatological 
averages of the pressures and adjusting by an S 2 model, 
here based on Ray and Ponte [2003]. 

[ 78 ] GOT4.9 is identical to 4.8, except that (unlike all 
previous GOT versions) it applied a correction to the T/P 
sea-surface heights to account for movement of the T/P 
center of mass relative to the altimeter antenna. The correc- 
tion is based on the orientation of the satellite, which 
changed primarily according to a yaw steering procedure 
with the sun/orbit-plane angle, as well as on the heating of 
the solar panel. (In the AVISO User Handbook this is the 
“CG_Range_Corr” correction, based on algorithm sl038; 
see also Callahan [2010].) The correction is controversial 
because of some anomalous effects seen during the Jason- 1 
calibration campaign when T/P and Jason- 1 were flying in 
tight tandem formation. 

[ 79 ] GOT4.10 is identical to 4.7, 4.8, and 4.9, except that 
only data from Jason- 1 and Jason-2 were used, with no T/P 
data. Otherwise, all other processing and data were kept 
unchanged. Note that the dry- tropospheric correction for 
the Jason satellites never had the S 2 air-tide error that the 
original T/P data had. 

[so] Finally, some of the GOT tidal solutions have been 
further adjusted after it was realized that altimetry can be 
usefully corrected for crustal loading deformation by the S 2 
air tide. The radial displacement, computed by standard 
methods [ Farrell , 1972] for loading from a global baromet- 
ric model [Ray and Ponte, 2003], is never more than 
1.4 mm, so the effect has been understandably neglected, 
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Table 6. Model Versions and RMS Differences With Tropical S 2 Bottom-Pressure Stations 



Used Topex 

Used Jason 

Corrected Dry-Trop 

Applied T/P C_g 

Corrected Air-Tide Load 

RMS (cm) 

GOT4.7 

Yes 

No 

No 

No 

No 

0.596 ± 0.040 

GOT4.8 

Yes 

No 

Yes 

No 

No 

0.371 ±0.023 

GOT4.8a 

Yes 

No 

Yes 

No 

Yes 

0.339 ±0.023 

GOT4.9 

Yes 

No 

Yes 

Yes 

No 

0.507 ±0.028 

GOT4.10 

No 

Yes 



No 

0.306 ±0.023 

GOT4.10a 

No 

Yes 



Yes 

0.273 ±0.021 


but it can influence results at the precisions now reached. 
The correction was applied only to GOT4.8 and 4.10, and 
these variants are here distinguished by an additional suffix 
“a” in their names. 

[8 1 ] Table 6 summarizes the major differences between 
these GOT solutions and also tabulates the rms differences 
between the S 2 solutions and the set of bottom-pressure 
tide stations described in the main part of the paper. Similar 
to Table 4, only tropical stations have been used because 
these effects are most pronounced in low latitudes. (The 
main conclusions apply also if the full set of 151 stations is 
used, but the rms differences are slightly more similar 
among models.) The tabulated standard errors were com- 
puted as above via a bootstrap resampling. 

[ 82 ] Several conclusions are evident in the rms statistics 
of Table 6. The altimeter-minus-pressure differences are 
clearly precise enough to show: (1) the revised dry- 
troposphere correction for T/P is very beneficial even 
though the original air-tide error was no larger than 
±3 mm; (2) the correction of satellite altimetry for crustal 
loading by the S 2 air tide is beneficial; and (3) the center- 
of-mass correction for T/P appears problematic. Less 
clearly, but probably significant because the effect is at the 
2(7 level, the differences between GOT4.8 and 4.10 again 
highlight an inherent inconsistency between T/P and the 
two Jason satellites, with the rms statistics perhaps suggest- 
ing Jason the more accurate. One might speculate whether 
another correction for T/P, similar but different from the 
center-of-mass correction now available, might bring T/P 
and Jason into better consistency, but it is unclear at this 
point how this might be done. The results do suggest that 
further studies of the T/P and Jason- 1 tandem campaign are 
still warranted. The inconsistency clearly affects present 
and future tidal modeling, since it is unclear how the two 
satellites should be weighted in altimetric assimilation and 
inversions. But it also contributes to the 59 day oscillations 
seen in time series of global mean sea level [e.g., Masters 
et al., 2012], since no one tide model — even one that 
absorbs other altimetric errors — can consistently remove 
the 59 day signals throughout the whole combined time 
series. 

Appendix B: Modeling Solid Tides in Altimetry 


[83] For the analyses in section 5, it is worth noting how 
the solid tide has been handled in the altimetric tide solu- 
tions. The solid tide is here considered comprising two 
parts : a body tide forced by the astronomical tidal poten- 
tial, and a smaller, more complex, load tide forced by 


crustal loading from the weight of the overlying ocean 
tide. 

[84] In nearly all altimeter tidal analyses the body tide is 
removed as a prior correction to the altimetric sea-surface 
heights. On both T/P and Jason Geophysical Data Records 
(GDRs) the body tide has been computed from a harmonic 
expansion of the tidal potential [ Cartwright and Edden, 
1973], extrapolated from 1960 to the 2000 era, with 
adopted real Love numbers h 2 = 0.609 and hi = 0.291, 
held fixed for all frequencies of degree-2 and degree-3 
tides, respectively, except that h 2 was reset to 0.52 for Kj 
and its two largest nodal lines. The latter reflects depend- 
ence of the h 2 Love number on the resonance effect of the 
fluid core [Wahr, 1981]; the GDR algorithm neglects even 
larger resonance effects at other constituents (e.g., ip i), but 
the body-tide displacement errors are small because the 
constituents themselves are so small. 

[85] The GDR Love numbers are not consistent with IERS 
Conventions [Petit and Lnzum, 2010]. The value 0.609 corre- 
sponds to results obtained by Wahr [1981] for semidiurnal 
frequencies, and the value 0.52 corresponds to his result for 
the main Ki line, but not the nodal lines nor any other diurnal 
line. Long time series of synthetic body tides, generated at 
latitude 30°N with the GDR algorithm and with the IERS 
algorithm, give an rms difference of 1 .2 mm ; the most diver- 
gent frequency is at Pi, where the core resonance effect is 
still significant (Wahr gives h 2 = 0.581 for Pi). 

[86] Unlike the body-tide correction, the load-tide cor- 
rection is a posterior one, determined after tidal analysis of 
the altimetry has yielded a combined ocean plus load tide. 
The load tide consistent with this solution can then be com- 
puted by an iterative method ( Cartwright and Ray, 1991, 
Appendix). For this, the loading numbers h 2 from Farrell 
[1972] were adopted. 
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